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Abstract 


The intensity and polarization characteristics of light scattered by bubbly 
media is studied, assuming that the size of bubbles is much larger than the 
wavelength of the incident light, which is a valid assumption, e.g., for 
bubbles in water and sea ice. The influence of possible nonsphericity of 
bubbles, light absorption and multiple light scattering effects is neglected. 
The dependence of the phase matrix on the relative refractive index of 
bubbles n is studied in the range n = 0.3-0.9. Results presented can be used 
for the development of analytical models of polarized light propagation in 
bubbly media, which take into account special features of bubble light 


scattering. 
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1. Introduction 


Optical properties of bubbles are of interest for many 
applications, including light propagation in oceanic water, 
ice and various technical media [1-4]. Usually the size 
of bubbles a is much larger than the wavelength of visible 
light A. Then phase function and polarization characteristics 
are independent of their size everywhere except the forward, 
backward, and critical scattering regions, where the light 
diffraction should be taken into account. For spherical bubbles, 
this can be achieved using the Mie theory [5], which correctly 
accounts for diffraction effects. Physical optics and complex 
angular momentum approximations [6-9] can be of help to 
see the physics behind a diffraction optics phenomena, which 
otherwise is hidden in the complicated Mie solution. 

The main aim of this paper is to study the dependence of 
the phase matrix of spherical bubbles on their refractive index, 
assuming that a >> A. The correct way to do this is to use the 
Mie theory for a given bubble size distribution [5]. We present 
such results for the case of air bubbles in water, assuming that 
the effective radius of droplets is equal to 15 and 30 um. 

It is known that geometrical optics can be considered as 
a limiting case of wave optics as à — 0. Therefore, it is of 
importance from a general point of view to derive these limiting 
patterns for elements of the Mueller matrix of bubbly media 
as à — 0 at various values of the refractive index n. The 
numerical calculation of these limiting patterns constitutes a 


major part of this work. This is essentially the extension of 
work [10], which was performed in the assumption that the 
relative refractive index of bubbles is 0.75. Only the intensity 
of scattered light and not a full Mueller matrix was considered 
in [10]. 

The modern ray tracing code, used in calculations, is 
described elsewhere [11]. The physical idea behind the 
geometrical optics approach, which is the base of code [11], is 
quite simple: the tracing of the geometrical optics rays as they 
interact with a bubble surface. The Fresnel equations [10, 11] 
are used to find the light energy and polarization change at 
the bubble—host medium interface. Results obtained with 
code [11] for the intensity of scattered light coincide with those 
given in [10] at a particular case of n = 0.75. 


2. Intensity and polarization characteristics of 
scattered light 


The interaction of light with a single bubble can be described 
in the framework of the Stokes vector formalism [12]. Then 
the Stokes vector of scattered light S is related to the Stokes 
vector of the incident light Sp by a following formula in a linear 
optics approximation [12]: 


qa 7 PO So, (1) 


where @ is the scattering angle, k = 27/2, r is the distance to 
the observation point and F (0) is the dimensionless scattering 
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Figure 1. The phase function of bubbles for various n = 0.3, 0.4, 
0.5, 0.6, 0.75, 0.8, 0.9, obtained in the framework of the geometrical 
optics. 


function, which defines the scattering cross section as [12] 


_ 2m 


Csca = K 


[ F(6) sin 6 dd. (2) 
0 


The normalized scattering matrix p, which coincides with the 
normalized Mueller matrix of the problem at hand, has the 
following general structure for spheres [12]: 


1 pm O 0 


à p2 1 0 0 
= . 3 
á 0 0 P44 P34 6) 
0 0 —p34 pu 
We note that 
27 y . 
C.K Jy F(@)sin0 dé = 1. (4) 


In what follows we study not the function F(@), but the so- 
called phase function 


An F(0) 


6) = —, 
p(@) CZE 


(5) 
which gives a conditional probability of light scattering in a 
particular direction 0 and normalized as follows: 


TOE (6) 
0 


Due to large size of droplets, the value of Cy. is equal to 
2S [2, 6, 12], where S is the average geometrical cross section 
of bubbles. 

We underline that equation (1) is also valid for 
spherical polydispersions. Then the function F(@), products 
F(0)pij (i, j = 1,2,3,4) and scattering cross section Csca 
should be averaged with respect to the bubble size distribution 
function f(a) [12]. 

The phase function p(@), calculated using the Monte Carlo 
ray tracing code [11], is given in figure 1 for various relative 
indices n of monodispersed bubbles (n = 0.75 corresponds 
to the case of air bubbles in water). The geometrical 
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Figure 2. The dependence of the asymmetry parameter g on n, 
obtained in the geometrical optics approximation (circles). Curve 
gives the Gauss fit results (see equation (8)). 


optics calculations were supplemented by the wave optics 
results at small scattering angles, using the well-known Airy 
formula [12] (F (0) ~ 4J? (ka0)/0° as 0 — 0, where J; 
is the Bessel function). We assumed in calculations that the 
radius of a bubble a = 1 mm and à = 0.55 um. So we 
have § = $ œ% 1818, which explains a sharp peak in the 
forward direction. Its width is proportional to A/a. As was 
stated above, p(@) does not change with £ (except at small 
scattering angles, if the diffraction is taken into account) in the 
framework of the geometrical optics approach. Thus, results 
given in figure | correspond to the asymptotic scattering regime 
(E > oo). 

The behaviour of phase functions is completely different 
for angles 0 < 0. and @ > 0., where 0. = 2 arccos(n) is the 
so-called critical scattering angle [2, 9]. The value of 6, grows 
with n — 1. Formally, we have 6. = 0° atn = 1 and 6, = 180° 
atn = 0. It follows at n = 0.75, 0, ~ 82.8°. 

In particular, we see that for n = 0.3 the phase function is 
almost constant in the region of scattering angles smaller than 
approximately 140° (with the exception of the diffraction zone 
at small angles). It sharply decreases at 0 = 0. ~ 145°. 
The phase function is again almost constant at 0 > 145° 
(except as 6 — 180°, which is due to the axial focusing 
effect [5]). Similar behaviour is observed for other n (see 
figure 1). Generally, the sharpness of the step at 0 = 6, 
increases with An = 1—n. The phase function p(@) coincides 
with the delta function ô(0) at n = 1. It is given by a 
combination of the Airy diffraction pattern at small angles and 
an isotropic part outside of the diffraction peak at n = 0. This 
is due to the fact that waves do not penetrate inside a particle 
in this case. 

The forward—backward asymmetry of phase functions is 
usually represented by the asymmetry parameter [12] 


CO 
g= if p(0) sin 6 cos 6 dé. (7) 
0 


This parameter is given in table 1 and figure 2 for various n. 
We see that g increases as n — 1. The solid curve gives a fit, 
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Figure 3. As in figure 1, but for pi2 (a), pas (b), p34 (c). 


using the Gaussian distribution 


ae A {—2("<") | (8) 
§ = 8&0 AA P A ’ 


where go = 0.492 74, no = 1.048 82, A = 0.692 33, and A = 
0.444 66. Note that the parameter g determines reflectancer = 
1 — tf and transmittance ¢ of optically thick bubbly layers [12]. 
In particular, we have for the large optical thickness t of a 
bubbly medium [12] 


1 


l= 10740750 — 2)" ©) 


Bubbly media with the same values of t(1 — g) have close 
values of r and t as T — ov, despite having quite different 
single scattering laws (see figure 1). 

The elements pj2, p44, p34 Of the matrix p, given by 
equation (3), are presented in figures 3(a)-(c). The element 
Piz (see figure 3(a)) is equal (with opposite sign) to the degree 
of linear polarization p; of scattered light, assuming that the 
incident light is unpolarized. Therefore, we see that p; is 
mostly positive (the oscillations of the electric vector are in the 
plane perpendicular to the scattering plane, containing incident 
and scattered beams). The degree of linear polarization sharply 
increases at the critical scattering angle. Then it takes a 


Table 1. The asymmetry parameter g for various n. 


n g 
0.0 0.5 

0.3 0.5380 
0.4 0.5807 
0.6 0.7160 
0.75 0.8443 
0.8 0.8869 
0.9 0.9614 
1.0 1.0 


maximum and again decreases for larger scattering angles. The 
maximum is sharper for smaller n. Note that there is also an 
increase in the degree of linear polarization in the vicinity of the 
backscattering direction. Also, we have p;(0) = p;(z) = 0 
for all n. This suggests that there is no light depolarization in 
the forward and backward scattering directions for spherical 
bubbles. 

The element p44, presented in figure 3(b), gives the degree 
of circular polarization pe of scattered light, assuming that 
incident light is right-hand circularly polarized (pe = 1). In 
particular, we see that pe = —1 for scattered light at exactly 
the backward direction for arbitrary n. This is due to the 
conservation of rotation sense (photon angular momentum) 
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in the light backscattering by bubbles (and spheres in general). 
Generally, light scattering by bubbles reduces the circular 
polarization of incident circularly polarized light (|p.| < 1). 
In particular, the value of p, is equal to zero for four scattering 
angles at n = 0.3 (see figure 3(b)). The strong variation of 
Dc is observed near the backward direction. For air bubbles in 
water (n = 0.75) pe varies very little with 6 < 6.. However, 
it drops sharply at the critical angle, reversing the sign of pe 
around 0 ~ 105°. Therefore, the scattering process at angles 
larger than 105° for air bubbles in water transforms incident 
right-hand polarized light to a left-hand circularly polarized 
mode. 

The element p34, presented in figure 3(c), gives the value 
of pe for the case of linearly polarized incident light with 
the azimuth —45°. We see that the value of p. is equal to 
zero at 0 > 0e. Therefore, linearly polarized light cannot be 
transformed in a circularly polarized mode at 9 > 6, in the 
geometrical optics scattering regime. 

Following this brief discussion of various matrix elements, 
we would like to underline that data given in figures 1 and 3 
(in combination with equation (1)) can be used to study the 
interaction of arbitrarily polarized light beams with bubbly 
media in the framework of the single scattering approximation, 
assuming that A/a — 0 [6]. 

It is interesting to see how the geometrical optics results 
described above correspond to the wave optics results. For this, 
we have presented the comparison of geometrical optics [11] 
and Mie [5] calculations at n = 0.75, A = 0.55 um in 
figures 4(a)-(d). Geometrical optics curves in figures 4(a)— 
(d) (solid curves) coincide with those given in figures | and 3 
atn = 0.75. The bubble size distribution in Mie calculations 
was assumed to be 


f(a) = Aa® exp(—9a/aer), (10) 
where a is the radius of droplets and A = constant. The 
function f(a) is normalized as follows [12]: 

oo 
| f(a) da = 1. (11) 
0 


We have used two values of the effective radius aef, namely 
15 and 30 um, in Mie calculations. 

It follows from figure 4(a) that phase functions calculated 
using geometrical optics and wave optics differ considerably 
around 0 ~ @, (and in the diffraction zone, which is mostly 
due to different size of bubbles used in wave and geometrical 
optics calculations). They depend on the size of particles 
and wavelength in the region of critical scattering [6-9]. 
This dependence is totally absent in the framework of the 
geometrical optics description [9]. We see that the scattering 
angle corresponding to the maximum of scattering in the 
critical scattering region increases with per = 27ae7/d. Its 
position approaches the asymptotic value as pef — oo. The 
wavelength dependence produces coloured rings in scattered 
light around 6 ~ @, if one uses white light illumination. 
Red light with larger A (smaller pef) gives a maximum at 
smaller 0 as compared with blue light for the same bubble 
size distribution (see figure 4(a)). The intensity of red rings is 
expected to be larger than blue ones. On the other hand, the 
critical scattering of a monochromatic light (e.g., laser) can be 
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used to estimate the value of aef [7]. Note that coloured rings 
in white light scattering patterns for bubbly media have been 
observed experimentally [13, 14]. 

It follows from figure 4(a) that the geometrical optics has 
a high accuracy in the backward hemisphere (9 > 90°) for air 
bubbles in water at pef larger than 180. The corresponding 
limiting value of pef for water droplets in air is larger (around 
1000 [15]). This is due to the rainbow effect [15], which cannot 
be described by a geometrical optics alone. 

The rainbow effect does not exist for bubbles, however. 
Therefore, main errors of the geometrical optics for small 
bubbles are in the vicinity of the critical (0 ~ 82.8°) and glory 
(0 ~ 180°) scattering regions, where the diffraction of light 
should be properly taken into account. The discrepancy around 
small scattering angle (9 ~ 0°) is mostly due to different values 
of the bubble radii assumed in the geometrical optics and Mie 
calculations. 

A close examination of figures 4(b)-(d) shows generally 
the same features as we found in figure 4(a). Namely, major 
differences arise in the vicinity of the critical and glory [6-9] 
scattering regions. The dependence of p12, p34, and p44 ON aef 
can be used for bubble optical sizing. 

It is interesting that the position and the value of the 
maximal polarization (see figure 4(b)) is described by the 
geometrical optics method with a high accuracy. This feature 
can be applied to measurements of the refractive index of a 
host medium with spherical air inclusions. 


3. Conclusion 


We have studied here the dependence of the normalized matrix 
elements pj2, p34, p44 and the phase function p(@) on the 
relative refractive index of bubbles with sizes much larger than 
the wavelength of incident light. The possible light absorption, 
multiple light scattering and change of frequency during the 
scattering process is neglected. We also ignored the possible 
distortion of droplet shape [16, 17] and the presence thin 
organic films, which can be formed on the surface of bubbles in 
many cases [14, 18]. These assumptions allowed us to derive 
an idealized geometrical optics pattern picture, which can be 
further clarified, if needed, to account for additional factors 
(e.g., using Mie theory). 

Note that bubbles can be air trapped in various materials. 
The medium of particular importance is, of course, ocean. In 
particular, oceanic whitecaps play an important role in the 
air—ocean exchange processes, modification of the planetary 
albedo and satellite remote sensing (as an obstacle for 
an accurate determination of aerosol optical thickness and 
chlorophyll concentration from orbiting satellites) [19]. It 
should be stressed that results obtained are generally applicable 
to any spherical particle having refractive index smaller than 
that of a surrounding medium (e.g. water droplets or bio- 
particles trapped in glass or various solids). 

The matrix (3) plays a fundamental role in the vector 
radiative transfer in multiple light scattering media [12, 20]. It 
can be used, therefore, to facilitate the derivation of analytical 
results in the field of polarized multiple light scattering, taking 
into account the peculiarities of single light scattering by 
particles with n < 1, considered here. 


10000 , 


100 | 
Cc 
S 
v ; 
§ w| 
oO 
i?) 
oO 
<= 
a 
i 
0.1 
0.01 


O 20 40 60 80 100 120 140 160 180 


scattering angle, degrees 


“0 2 40 60 80 


100 120 140 160 180 
scattering angle, degrees 


Optical properties of bubbles 


b) 


Piz 
<4 
ro) 
| 


O 2 40 6 80 100 120 140 160 180 
scattering angle, degrees 


0 20 40 60 80 100 120 140 160 180 


scattering angle, degrees 


Figure 4. (a) As in figure 1, but at n = 0.75 only (solid curves). Dashed curves give the Mie theory results for air bubbles in water, having 
the size distribution (10) with a.p equal to 15 zm (long-dashed curve) and 30 um (short-dashed curve). The wavelength is equal to 
0.55 um. (b) As in (a) except for pı2. (c) As in (a) except for p44. (d) As in (a) except for p34. 


We underline that the critical scattering region can be used 
for optical bubble sizing [7]. This could be a useful addition 
to small-angle particle sizing methods [12]. 
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